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1. Introduction 



The presence of a gluon and quark 'condensate' in the QCD ground state is a plausible 
reason for the observed long distance properties of QCD p], 0. The condensate appar- 
ently prevents quarks and gluons from propagating over long distances, while acting as a 
superfluid for color singlet hadrons. In this work we model the gluon condensate effects 
by coupling quarks and gluons to a 'vacuum' gluon field <3? which is taken as a constant 
in space-time in a covariant gauge. Translation invariance is thus automatically satisfied. 
We integrate over the Lorentz and color components of ^ with a gaussian weight. This 
ensures Lorentz and gauge invariance and introduces a dimensionful parameter A which 
characterizes the magnitude of the vacuum field. In effect, we modify QCD perturbation 
theory (PQCD) by expanding around non-vanishing gluon field configurations. 

In the limit of a large number of colors N —* oo with g 2 N fixed Q, we are able to find 
the exact expressions for the gluon and (massless) quark propagators in a 'dressed tree' 
approximation. The vacuum field effects are resummed to all orders, whereas perturbative 
loop corrections are neglected. The dressed tree propagators have a branch cut 

instead of a pole at p 2 = 0, and consequently decay with time t as l/yt. For \p 2 \ 3> 
fj? = g 2 NA 2 the dressed propagators approach the free ones. Hence the short distance 
structure of PQCD is unaffected by the vacuum field. It is gratifying that the momentum 
dependence of our dressed quark and gluon propagators in the large N limit turns out to 
agree qualitatively with the results of lattice calculations at iV = 3. 

The finite propagation length of the dressed partons appears to regularize the infrared 
(IR) singularities of the standard PQCD expansion. We study the dressed (massless) quark 
loop contribution to the photon self-energy. Zero-momentum gluons can couple to the quark 
loop even for spacelike photons due to the appearance of infrared singularities in Feynman 
diagrams involving the coupling of at least 4 external fields. While the contribution of a 
specific number of external fields is thus ill-defined, the loop integral is IR regular when 
the fully dressed quark propagator and quark-photon vertex are used. The effective lower 
cut-off for the loop momentum k is k 2 ~ /x 2 = g 2 NA 2 (in euclidean space). Thus the 
dressing indeed 'confines' the color singlet quark pair within a distance ~ l//i. At high 
photon virtualities p 2 the dressing correction is oc l/p 4 as in the QCD sum rule framework 
Q, where the normalization would be given by the vacuum expectation value (a s F^ u Fa U ). 

Since the elementary constituents of QCD are confined their Green functions need not 
have a standard analytic structure. Our dressed quark and gluon propagators indeed have 
branch cuts in p 2 . The present framework may thus give insights into how an S-matrix 
consistent with general principles can be constructed in a confining theory. A first step in 
this direction will be to identify the asymptotic states of our framework. 

This work represents a further development of our previous work where some of 
the results on the dressed quark propagator and photon self-energy were already presented. 

2. Coupling quarks and gluons to the vacuum field 

Gauge invariant couplings of quarks and gluons to the vacuum field $ can be found by 
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shifting the gluon field, 

A^-tAP + V (2.1) 

The quark part tpilftip of the massless QCD lagrangian (with = d^ + igA^ the covariant 
derivative) then generates the coupling 

U q = -g^i^ (2.2) 

which we shall use. It is invariant under the gauge transformation tp — * U{x)ip with 
U(x) £ SU(iV) provided $ transforms as 

^ -» U{x)^U{x) ] (2.3) 

This transformation law is consistent with that of the shifted gluon field, 

A» + ^ -» U{A^ + Q^rf + - (d"U) U ] (2.4) 



Under the shift (2.1) the field strength tensor 

F^ = d^A u -d v A tl + ig[A^A v ] (2.5) 

becomes 

Fnv -» + F%, + ^ (2.6) 

F% = d^ v - + ig{[%, A v ] - [*„, Afj]) (2.7) 

= ig *„] (2.8) 

since -> l/i^U* under a gauge transformation, each term in the r.h.s. of (|2.6j ) of a 
given degree in <J>^ transforms similarly: 

F^^UF^rf ; F^^UF^ ; ^ _» L^C/t (2 . 9 ) 

We wish to study the influence of a low-momentum 'vacuum' gluon field on quark and 
gluon propagation without delving into the more difficult question of how the condensate 
field is generated. We therefore ignore the self-interactions of <1> and consider only the 
gauge invariant coupling of the <3? field to gluons which is linear in 

£* g = -Tr [F^F%] (2.10) 

We shall thus work with the modified (massless) QCD lagrangian, 

C = $i$i/>- ^Tr [i^FH - |Tr [{d^f] - cd^D^c -Xg^^-Tr [F^F%] (2.11) 

which includes covariant gauge fixing and ghost terms and is BRST invariant. Since the 
quark and gluon couplings to $ are separately gauge invariant there is no constraint on 
their relative weight A. For simplicity we choose A = 1 in the following. 
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As the field <& is meant to describe the long wavelength (vacuum) effects we take it 
to carry zero momentum, i.e., to be independent of the coordinate x (in the covariant 
gauge specified by ( |2.1lD ). Translation invariance is thus guaranteed. In order to preserve 
Lorentz (and gauge) invariance we average over all Lorentz (and color) components of 
with a gaussian weight 1 , 



n 



d$£ exp 



1 



2A 2 " 6 



(2.12) 



where A is a parameter with the dimension of mass. In a perturbative expansion we may 
interpret ( |2.12| ) as giving a <I> 'propagator' 



iD 



ab 

<5>,fiu 



(2.13) 



In the next section we derive explicit expressions for the quark and gluon propagators 
in a 'dressed tree' approximation, which takes all interactions with the $ field into account, 
but neglects perturbative quark and gluon loops. 



3. Dressed quark and gluon propagators in the large TV limit 



We now calculate the effects of the interaction terms (2.2) and ( [2.10D of the zero-momentum 



vacuum field on the quark and gluon propagators. We simplify the topology of the 
contributing diagrams by taking the limit of a large number of colors, N — > oo with g 2 N 
fixed Q. In addition to the coupling g 2 N we have a parameter \x with the dimension of 
mass, 

f i 2 = g 2 NA 2 (3.1) 



With the lagrangian ( 2.11| ) the full perturbative expansion of any quark and gluon Green 



function G is a double sum of the form 

oo oo 

G = £GW£<W fc (3.2) 

1=0 k=0 

Here £ counts the number of perturbative loops and 2k the number of f> couplings ( |2.2j ) 



or ( 2.10 ). We evaluate the complete sum over k for £ = 0. This is possible since the <3? 
propagator ( 2.13j ) carries zero momentum so no loop integrals are involved. The remaining 
sum over £ involves higher powers of the coupling g 2 N and an increasing number of non- 
trivial loop integrals 



2 



3.1 Dressed quark propagator 

Only planar $ 'loop' corrections to the free quark propagator contribute in the large N 
limit (with £ = in ( |3.2| )). Prom the structure of a general diagram it is readily seen 



1 The integrals over the time components Of} are defined by analytic continuation. 

2 We do not consider quark and gluon loops which are connected to the rest of the diagram only via $ 
lines. These would generate self-couplings of the $ field. 
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that the dressed quark propagator S(p) satisfies the Dyson-Schwinger (DS) type equation 
shown in Fig. [l], which reads 



iS(p) 



+ C F (-ig) 2 (-A 2 )- r u iS(p)'y ll iS(p) 



(3.3) 



where /i is defined in (3.1) and we used Cf = N/2 at leading order in N . The DS equation 
( |3.3| ) generates all one-particle irreducible as well as reducible planar diagrams. 
Lorentz invariance constrains the quark propagator to be of the form 



S(p) = a p tf + b p 

where a p and b p are functions of p 2 . Substituting this in ( |3.3D we get 
6(1 + //a) = and ap 2 {\ - fi 2 a) = 1 - 2\i 2 b 2 



(3.4) 



(3.5) 



A chiral symmetry conserving quark propagator S\(p) must have 6 = 0. The second order 
equation for a then gives 



p 2 



P 2 + \Jp 2 {j) 2 — 4:{i 2 



(3.6) 



(3.7) 



where we chose the sign of the square root so as to ensure that the propagator approaches 
the free one in the p 2 — > oo limit: 



Si(p) 



l + O 



for p 2 



oo 



(3. 



The p 2 = pole of the free quark propagator i> jp 2 has been removed by the dressing. The 
dressed propagator S\(p) given in ( |3.7|) has instead branch point singularities at p 2 = 
and p 2 = 4/i 2 . In appendix A (see (|A.2|) ) we show that this removes the quark from the 
set of asymptotic states, in the sense that the Fourier transformed propagator vanishes at 
large times, 

(3.9) 



l^p)! ~ 0(l/y/\t 

\t\— »oo 



S{p) = 



P 



Figure 1: Implicit equation for the dressed quark propagator S(p). The dashed line denotes the 
$ 'propagator' (2.13). 
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Figure 2: Numerical results for the quark propagator in the euclidean region p 2 E = —p 2 > 0. The 
quantities M and Z are defined in (3.12). 

Allowing chiral symmetry breaking (XSB), i.e., b ^ in ( |3.5| ) we find 

S 2 (P) = ~ (j ± VP 2 + M 2 /2) = — ,\ ■ = (3.10) 

This solution is singular for fi 2 — ► and hence does not have a power expansion in [i 2 of 
the form ( |3.2| ), It emerges as a 'non-perturbative' solution of the implicit equation (p.3j). 
Like the chiral symmetry conserving solution Si, the propagator 52 has no quark pole, 
only a branch point at p 2 = —fj, 2 /2, where it coincides with Si. Since the solution 52 (p) 
does not approach the perturbative propagator \jf at large p 2 it must be discarded at 
short distance. 

In Euclidean space we can estimate the chiral symmetry breaking effect of using the 
solution S2(pe) for < p\ = —p 2 < /i 2 /2 and the S± propagator for p 2 E > fi 2 /2. The value 
of the quark condensate obtained in this way (using the upper sign in ( |3. 1Q| ) ) is 

m =Jj^Tr [S 2 ( PE )] 6(^/2 - pi) = (3-11) 

This estimate, \ipip) — —[i/9A, is an order of magnitude below the generic scale 
Expressing the quark propagator as 

sip) = j^m (3 - 12) 

we plot in Fig. [2| the functions M(p 2 ) and Z(p 2 ), using \i = 500 MeV and the above 
combination of S\ and 52 in the euclidean p 2 < region. Their shapes are similar to the 
results obtained in lattice calculations, see for instance Refs. pj or Fig. 3 of Ref. ||. 

Let us mention that the above solutions 5i and 52 for the quark propagator were 
previously obtained in a different framework Q , and recently within a model assuming the 
presence of an (A®A%) condensate 0]. In the rest of this paper we shall use S(p) = S\(p) 
for all values of p 2 . 
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3.2 Dressed gluon propagator 

It is convenient to use the double line color notation [g] for the fields and <j?^ in the limit 
N — » oo. The doubly indexed fields are defined in terms of the color generator matrices 
T a of the fundamental representation as 



#3 



arpa 3 



(3.13) 



In the covariant gauge defined by (2.11) the free propagators and their graphical represen- 
tations are then 



i I 

; k 



P =o 



i 

k 



2 k p 



(3.14) 



(3.15) 



The interaction term ( 2.1C| ) couples $ to two and three gluons, via &AA and $j4j4j4 
vertices, respectively. The <3?j4j4j4 vertex does not, however, contribute in the dressed tree 
approximation (the I = term in (|3.2| )), since it implies at least one perturbative loop 
integral. The $^4^4 vertex is 



Uaa = ~2igTr ( [%,A V ] (&>A V - VA»)) 

= -2igTr [A v , d»A v \ - $ M [A", d v A v \ ) 



(3.16) 



where we used dpQ v = and dropped a total derivative to obtain the second line. The 



second term of ( 3.16 ) vanishes in Landau gauge, d v A v = 0. However, we keep both terms 
in order to see the dependence on the gauge parameter £. In the double line notation the 
corresponding Feynman rules are: 



A } 



k 'W 1 n 



Ap n 

t 3 m 



(3.17) 




+4igp P gpuS^ 5™ 5^ 



(3.18) 



The rules for vertices where the $ field attaches to the lower gluon color line are the 
same (up to an overall sign). Since nonplanar diagrams like the one shown in Fig. ||a do 
not contribute at large N, the $ lines never cross the gluon propagator. A generic diagram 
contributing to our dressed tree approximation is shown in Fig. |3|b. 



-6- 



(a) 



(b) 



Figure 3: Diagrams contributing to the gluon propagator in the double line color notation: a 
non-planar subleading contribution (a) and a generic dominant contribution (b) in the large N 
limit. 



One might expect that the dressed tree gluon propagator would satisfy a DS equation 
(or a finite set of equations) analogous to the one satisfied by the dressed quark propagator 
(Fig. 1). However, because planar corrections to the gluon propagator can appear on both 
sides of the gluon line the dressed §gg vertex is not proportional to the bare vertex, as is 
the case for the quark. This apparently precludes a finite set of closed DS equations for the 
gluon propagator (which in our approach would imply an algebraic solution) . By explicitly 
summing all diagrams we indeed find the non-algebraic expression ( |3.23j ) below. 

We first do the calculation in Landau gauge (£ = 0), and then evaluate the £ depen- 
dence. In Landau gauge the free gluon propagator is transverse, 



^£>V i = k' V 4 Pf(p) ; Pfip) =r-^- (3.19) 

z p p 



and thus the circled vertex ( |3.18| ) does not contribute. Independently of the topology of a 
diagram, adding one (double) line always yields the same factor, namely, 



A 2 W-n 2 2fi 



^(-4* 5 )V^— =-^ = * (3-20) 

Let n(k) be the number of distinct planar graphs having k & lines (cf. (|3.2;)). The dressed 
gluon propagator in Landau gauge is then (with standard gluon color indices): 



iD ab(p) = ^ P f(.P)d(-^-)Saj, with d{x) = y^n{k)x k . (3.21) 



k=0 



In appendix B we show that n(k) = C k Cfc+i, where Ck is a Catalan number (see (|B 
and ( |B.5[ )). Mathematica then gives 3 for the 'gluon polarization function' d(x): 



d{x) = f^C k C k+1 x k = l ^ 16X) with h(x)= 2 F x (-\, \, 2,xV (3.22) 

b— n \ / 



k=0 

The dressed gluon propagator is thus, in Landau gauge, 



i / i 1 1 2 ^ 



lab 



(3.23) 



3 The expression (3.22) for d(x) may be verified by noting that h(x) satisfies the differential equation 
x(l - x)h" + (2 - x)h! + A/4 = with the condition h(0) = 1. 




Figure 4: Real part (solid line) and absolute value of the imaginary part (dashed line) of the gluon 
polarization function d{x) in the euclidean region p 2 E = —p 2 > 0. 

From the integral representation of the hyper geometric function h(x), 



h(x) = - du y/l - v? \J\-xu 2 (3.24) 
vr In 



we see that h(x) has a branch cut for x > 1. This implies that the gluon propagator ( 3.23) ) 
has a cut for 

-32^ 2 <p 2 <0, (3.25) 

which, surprisingly, lies in the spacelike p 2 < region. 

The polarization function d{x) approaches unity in the limit of high gluon virtuality, 
x = -2fi 2 /p 2 -> 0: 

d( x 0) = 1 + 2x + Wx 2 + O (x 3 ) (3.26) 

Thus the dressing effects are power suppressed at short distances, as expected. In the long 
distance limit p 2 — > (x — > oo) the function d(x) vanishes, 



which turns the p 2 = pole of the free propagator into a l/y^p 2 " singularity. As for the 



quark, this implies that the gluon propagator decays with time t as l/Vi, see ( A.3 ) of 
Appendix A. 

We plot d{—2y?/p 2 ) for p 2 < and fi = 500 MeV in Fig. ||. The real part has a shape 
similar to the gluon propagator found numerically in a Landau gauge lattice calculation, see 
Refs. 0] or Fig. 3 of Ref. ||. However, the presence in our expression ( [3.23 ) of a branch cut 
in the euclidean region p 2 < gives a negative imaginary part to the polarization function 
which obscures the comparison. 



Our result ( [3.22 ) for the polarization function d(—2fj, 2 /p 2 ) was derived in Landau gauge 



(^ = in ( |3.14D ). We now show that this expression actually is independent of £, i.e., the 
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dressed gluon propagator in a general covariant gauge is 



-% 
P 2 



(3.28) 



The ^-independence of d(x) actually holds separately for any diagram contributing to 
d(x), planar or non-planar. This can be seen as follows, using the standard single line color 
notation for the gluon and propagators. Let k be the number of <J>-lines of a diagram 
with given topology. Since the free gluon propagator has a longitudinal part oc £, 

«3&(p) = = j ( p f(p) + tP£T (p)) ^ ; pT(p) = ( 3 - 29 ) 

the circled vertex ( |3.18| ) can contribute when £ ^ 0. There are 2 2fc contributions for a given 
topology, since a chosen vertex can be either circled or not. When no vertex is circled the 
Lorentz structure brought by the gluon line is oc (Pt + £,PL) 2k+1 = Pr + £ 2?c+1 -Pl- When at 
least one vertex is circled, consider the first vertex of this type one meets by following the 
gluon line (carrying momentum p) from left to right. In all gluon propagators to the right 
of this vertex only the longitudinal part Pl contributes. The m gluon propagators to the 
left of this vertex give (-Pt + C-Pl)" 1 = Pt + ^Pl, *-e. these m propagators are of the same 
type, either transverse or longitudinal. Independently of the topological configuration, 
(Pt + £ m Pi J )t il '' is always contracted with p v i [y 1 being the Lorentz index of the gluon line 
just before the first circled vertex). As a result, in all (2k + 1) gluon propagators, only 
the longitudinal part contributes when at least one among the 2k vertices is circled. Since 
the circled vertex has a negative sign relative to the uncircled one, the ^-dependent part 
brought by the 2 2k contributions is proportional to: 

f +1 F L Q-l)"p fc ]=0 (3.30) 

n=0 V n / 

Only the ^-independent (transverse) contribution remains, which is obtained when no ver- 
tex is circled. Thus the result found in Landau gauge for k > 1 holds in all covariant 
gauges, establishing fl3,2§| ). 



4. Photon self-energy 

In the previous section we saw that the dressed quark and gluon propagators have no pole 
at p 2 = 0, only a type singularity. This implies a limited propagation length and 

should soften the collinear (p 2 — > 0) and infrared (p — > 0) singularities of the perturbative 
expansion. In this section we study the effects of dressing the quark loop contribution to 
the self-energy of the photon. 

4.1 Coupling of zero- momentum gluons to a color singlet quark loop 

Spacelike photon fluctuations j(p) — > qq — ► j{p) involve quark pairs of size ~ p 2 
(for massless quarks). One might expect that such compact, short-lived color singlet states 
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would decouple from gluons of zero momentum, i.e., from the field we introduced in 
section |2[ This charge coherence effect is most easily seen for QED amplitudes. The 
standard contribution of a bare electron loop to the photon self-energy is 



U a0 (p) 



d d k 
(2vr)° 



Tr 



^S e (k)j a S e {k) 



where k = k — p. It is instructive to retain the mass in the electron propagator 

1 



S e (k) 



m 



(4.1) 



(4.2) 



Now consider the attachment of an external zero-momentum photon with Lorentz index 
Hi to this loop. This contribution, denoted by IT^f (p), is given by two diagrams. Using 
the identity 

" ' ' ^ : ' (4.3) 



m 



m 



7w 



m 



we note that IT^f (p) is obtained from (4T) by differentiating the integrand with respect to 
k^\ 



<(P) 



ze 



f d d k 


d ' 






Tr 


1 (2ir) d 


e dk^_ 



jPSeik^Seik) 



(4.4) 



Similarly, the QED amplitude for an electron loop with two external zero-momentum pho- 
tons is 



n a/3 



(p) 



te 



d d k 
(2vr) d 







dk^ 



d 



dk^ 



Tr 



^S e (k) 7 a S e (k) 



(4.5) 



In Eqs. (|4.4|) and (4.5) the integrand is given by a total derivative. This generalizes to 
any number of external zero-momentum photons. Thus the loop integral formally vanishes, 
in agreement with intuition that zero-momentum photons do not couple to a virtual e + e~ 
pair in a photon. The analogous proof for QCD is somewhat more involved and is given in 
Appendix C. 

However, this demonstration fails if the integral is singular and thus ill-defined. As 
seen from ( |4.3| ) the insertion of a zero-momentum photon increases the number of electron 
propagators having the same momentum, which may cause infrared singularities. Consider 
the photon self-energy contribution with n zero-momentum photons attached to a (massive) 
electron loop. In this amplitude there are diagrams with up to n + 1 electron propagators 
of the same momentum, e.g., [l/(^ — m)] n+1 . After a Wick rotation the integral has 
contributions at low k of the form (for even n) 



d A k¥ 



k-k 



{k 2 + m 



2^71+1 



(4.6) 



For a finite electron mass m ^ the integrand is always regular at k = and, being a 
total derivative, the loop integral must indeed vanish. On the other hand, for m = the 
integral is apparently IR regular only for n < 2 (and thus vanishes for n = 2 when the 
UV behaviour is dimensionally regularized). For n > 4 the loop integral is IR singular and 
hence ill-defined. 
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r^(k,k) = [L P 





Figure 5: Implicit equation for the quark-photon vertex r A '(fc,fc). 

The physical reason for this infrared behaviour is that the coupling of a photon to 
an e + e~ pair is proportional to the electric dipole moment of the pair, which favors large 
dipole configurations in the integrand. For n > 4 the integral becomes divergent at k = 0. 
The electron is then delocalized in space and by the uncertainty principle the pair can have 
any size. At k = p the positron is similarly delocalized. This long-distance behaviour is a 
physical effect which cannot (or should not) be removed, e.g., using dimensional regular- 
ization. In the case of a finite electron mass the maximum size of the virtual pair is set by 
1/m and no infrared singularities occur. 

This analysis also applies to QCD. Zero-momentum gluons can couple to virtual qq 
states in the photon and the dressing affects the photon self-energy. Using our dressed 
(chirally symmetric) quark propagator ( |3.7D , and the dressed 7qq vertex which we shall 
next derive, we find an infrared regular loop integral. The dressing thus keeps the quark 
pair at a finite separation ~ 1/fJ,. Singularities only appear if one tries to Taylor expand 
in powers of /i 2 , implying that non-analytic terms in fi 2 appear in a [i — > expansion. 

4.2 Dressed quark-photon vertex 

In our dressed tree approximation (I = in Q3.2| )) the 7qq vertex T ll (k,k) satisfies (at 
large N) the implicit equation of Fig. ||, 

r"(fc, jfe) = y* - ±n 2 YS(k)T»(k, k)S(k)j p (4.7) 

where k = k — p and p is the photon momentum. When the vertex is expanded on its 
independent Dirac components (|4.7l ) reduces to a set of linear equations for the compo- 



nents (c/. Eqs. ( D.4 ) and ( p.7| )), with a unique solution when the quark propagator S is 



given. The explicit expression (|D.9[) for T^(k, k) is derived in Appendix D using the chiral 
symmetry conserving quark propagator (|3.7|) . 



Multiplying (4/7) byp M = {k — k)^ and rewriting (3.3) as 



^VS^^^)- 1 -^ (4.8) 

it is readily verified that the solution of ( f4,7D respects the Ward-Takahashi identity 

P^(k,k) = Sik)' 1 - Sik)" 1 (4.9) 

For highly virtual momenta k 2 — > oo (at fixed p), we have S[k) —> S(k) — ► 1/k 
and ( |4.7p thus implies 

k 2 ^oo T^{k,k)=^ + o(^\ (4.10) 
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This may also be verified from the explicit expression ( p.S| ) of the dressed vertex. 



4.3 Dressed quark loop 



The dressed quark loop correction to the photon propagator is given by the dressed quark 
propagator and quark-photon vertex as indicated in Fig. 0. In terms of the solutions of 



the DS equations for the quark propagator ( |3.3| ) and vertex (4.7) the photon self-energy 
correction Ii pu (p) is 

d d k 



ie 2 N 



2 ie 2 N 



d-2 



where we used the relation 

k 
k 



Figure 6: Photon self-energy IF 1 "^) 
dressed by the <& field. The dressed prop- 
agators and vertex are indicated by a 
solid circle. 



d Tr [YS(k)T»(k,k)S(k)] 



d d k 



Tr[(r^(M)- r)l V \) 



(4.11) 
(4.12) 




li" J {2-K) d 

for the vertex and 7^ = j p r y u j p /(2 — d) in d dimensions. 

We dimensionally regularize the standard ul- 
traviolet divergence which appears at zeroth order 
in fi 2 . At O (/x 2 ) the quark loop couples to two 
photons and two $ fields and is UV finite, but is 
still ill-defined in d = 4 because of the superficial 
logarithmic divergence of the loop integral. Due to 
( |4.5| ) it vanishes when dimensionally regularized. 

As is well-known, within dimensional regular- 
ization the Ward-Takahashi identity (^1]) implies a 
transverse photon self-energy 



U^(p)=U(p 2 )(p 2 g 



2 fJIV 



p»p v ) 



and a resummed photon propagator 



P 2 [i-u( P 2 )} 



fJ.U 



pf i p u 
p2 



(4.13) 



(4.14) 



where p 2 H(p 2 ) — ► for p 2 — > 0. 

4.4 The dressed photon self-energy in euclidean space 

Here we consider in more detail the expression for U(p 2 ), specified by (gig ) and (JDJ), 
using the quark propagator S\(p) given in Q3.7| ) and the quark-photon vertex ( p.9| ). To 
avoid UV divergent terms in the loop integral we subtract the terms of O ux°) {i.e., the 
standard PQCD expression) and O (/i 2 ) (which vanishes in dimensional regularization) : 



U(p 2 ) = U(p 2 ) - U(p 2 



8ie 2 N f d A k 



3p 2 



n (p 2 )U 



H 2 k 2 k 2 ala 



(4.15) 



k ■ kakaj, + \i k k, ayx 



(2tt)' 



1 



^k 2 k 2 a\a\ 



1 + 2y?k ■ kakaj, + /z 4 /c 2 /c 2 a|a| 



+ 



fc-fc 

m 2 



p 
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where the function a p is given in (£ 

II(p 2 ) has a complicated analytic structure due to the denominators of the integrand 
and the square roots in a k and a^. We shall here only address its properties for euclidean 
momenta p 2 < 0, assuming that the loop integral can be Wick rotated. 

The denominator of the second term in the integrand of (4.15) can be estimated using 
\k-k\< \k\\k\, 



1 — 2fi 2 k ■ ka k a k + ^ 4 A; 2 /c 2 a|a| > (l — fJ- 2 \k\ \k\a k a k ) > 



(4.16) 



where we reversed the sign of all dot products since the momenta are now understood to 
be euclidean. This denominator can thus vanish only if 

4/272 2 2 

fx k k a k ar 



1 



(4.17) 



which is also the condition for the denominator of the first term in (4.15) to vanish. For 
euclidean k 2 > the function a k is negative, 



a k 




< 



The identity 



2 2 ; 2 

fi a k k 



1 + k a k 



(4.18) 



(4.19) 



then shows that fi a k k < 1 (with a similar relation for k). Consequently the condition 
(4.17) cannot be fulfilled for finite euclidean momenta, and the denominators in the ex- 
pression (4.15) for n are positive definite. 

It is remarkable that the loop integral in ( 4.15 ) is IR convergent. For k — > 0, we 
have k — ► —p and a k ~ —l/fiy/k 2 . Thus the first two terms of the integrand approach a 
constant, while the last two are IR safe. As we noted above in ([O]), the coefficients of \i n 
in the integrand become progressively more IR singular as n increases. The sum over all n 
implied by the dressing nevertheless gives a finite result. 

The O (/i 4 ) contribution to II(p 2 ) in ( [4.15 ) is 

16e 2 iV A 



n(^) 



3p 2 
d*k k 



(4.20) 



k 



[k A + k 4 -2k- k(k 2 + k 2 ) 



k 2 k 2 



2(k ■ kf 



(2tt) 4 k 6 k G 

The integral is UV convergent, but has superficial linear and logarithmic divergencies for 
k — ► (and similarly for k — ► 0). The linearly divergent terms are odd in k and thus do 
not contribute to the integral. The logarithmically divergent terms give 



U(p 2 



2- 



16e 2 N 



3p 2 



<Pk 1 
(2vr) 4 fcV 



1 



(fc • pf 
k 2 p 2 







(4.21) 



which vanishes due to the angular integration in d = 4. Hence the O (/U 4 ) term ( |4.20| ) is 
actually finite, and gives the leading behaviour of II(p 2 ) in the limit p 2 — > oo, 



n(/ 



7aN /j 4 
3tt p 4 



(4.22) 
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We recall that the defining expression ( [4,15 ) for II (p 2 ) is IR regular. A logarithmic 
singularity in the O (/i 4 ) expression ( f4.20| ) would have signalled an asymptotic behaviour 
oc (/i 4 /p 4 ) log(/x 4 /p 4 ). Since according to ( |4.6| ) the IR behaviour becomes more singular 
with the power n of \x it is likely that the next-to-leading term in the p 2 — > oo limit vanishes 
more slowly than O (^ 6 /p 6 ) • 



5. Concluding remarks 

We have formulated a modified perturbative expansion of QCD, where the standard di- 
agrams are dressed by a constant external gluon field <I>^ which is gaussian distributed 
in magnitude. We derived explicit expressions for the dressed massless quark and gluon 
propagators, as well as for the photon-quark vertex, at lowest order in the standard loop 
expansion. We discussed some general properties of the quark loop contribution to the 
photon self-energy. 

Our formulation is explicitly gauge and Poincare invariant, and at short distances 
introduces only power suppressed corrections to the standard PQCD results. The dressed 
propagators have a branch point singularity instead of a pole at p 2 = 0. Hence quarks and 
gluons cannot appear as asymptotic states, in accordance with intuition that colored objects 
do not propagate to infinite distances in a color field. We also found that the dressed quark 
loop contribution to the photon self-energy is regular in the infrared euclidean domain. This 
contrasts with the IR sensitivity of bare loops with four or more soft external gluons, and 
implies that the physical size of the dressed qq pair is governed by the scale 1/fi related 
to the strength of the field <I>. 

We introduced the external field $ as a means to study the qualitative effects of a 
gluon condensate on quark and gluon propagation. While the features mentioned above 
are encouraging, much work remains to be done to ascertain whether this method yields 
results which are in accordance with general principles. Our dressed quark and gluon 
propagators have a novel analytic structure. It is especially surprising that the dressed 
gluon propagator has a cut in the spacelike p 2 < region, and it will be crucial to check that 
this behaviour is not in contradiction with causality requirements. The analytic properties 
of Green functions for confined fields are largely unknown, an issue which the present 
framework may help to clarify. It will also be important to identify the asymptotic states 
(if any!) of our framework, and to check whether an analytic and unitary S-matrix can be 
defined. 
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APPENDIX 



A. Asymptotic behaviour of the dressed quark and gluon propagators 

We give in this appendix the asymptotic time behaviour of the dressed chirally symmetric 



quark propagator S±(p) ( |3.7| ) and of the gluon propagator in Landau gauge ( 3.23 ). 
The Fourier transformed propagators are defined by 



We find 



S(t,p) = [°° ^-S(p)e~ u ^ ; D>»(t,p) = f°° ^-D^(p)e~ Up ° (A.l 

exp(-i|tp|) 



S(t,p) 
D^p) 



1 + i 

r^i 

|t|— »oo 2y/7T 
2^2 



|p|7° - P-7) 



\t\— >oa T^\[t^ 



exp(— i|tp| 



(A.2) 
(A.3) 



Since the two calculations are similar we give only the derivation of ( |A.2| ) for the quark 
propagator. 



The quark propagator (3/7) may be written as 

p 2 - 4/x 2 



Si(p) 



2fx 2 



1 



a/p 2 + ie^/p 2 — 4^ 2 + ie 



(A-4) 



where the zs prescription arises from the usual Feynman prescription of the free quark 
propagator p/(p 2 + ie). 

The Fourier transform gives 



S(t,p) 
J{t,p 2 ,/j 2 ) 



{l7 ° dt 2 P ' f) [«(*) + (P 2 + 4^ 2 + d 2 ) J(t, p 2 ^ 2 ) 



2/x 2 
dp 



exp(-itp ) 



2tt 



\/Po - P 2 + ^\/Po - P 2 - 4 ^ 2 + ie 



The function J(t,p , /j, ) can be evaluated using Feynman parametrization, 



1 



1 



dx 



1 



\M + ie\AB + ie Wo y/^TTT"^) (1 - X )A + xB + ie 
and doing the j?cr m t e g ra l using Cauchy's theorem. The result is 



(A.5) 
(A.6) 

(A.7) 



J(t,p 2 ,^ 2 ) 



-i|tp 



dx 



2in Jo y/ x (l - x) a/p 2 + Ax^j 2 



cxp 



-4i|i|x/i 



a/p 2 + 4x[i 2 + |p| 



(Ai 



The behaviour of J(t,p 2 ,/x 2 ) for |i| — > 00 can be inferred by noticing that the integrand 
in ( |A.8| ) is peaked at x — > in this limit. With the change of variable y = 2\t\/j> 2 x/\p\ we 
obtain 



J(t,p 2 ,/x 2 ^ 



-i|tp 



dy 

zvr-^/Sltpl/x 2 Jo s/V 



-iy 



1 + i e — <[*■»! 



4aA v^pP 



(A.9) 



- 15 - 



where we used 



Using ( |OD in Q gives the asymptotic time behaviour (|A.2|), 



B. Counting planar graphs 

In this appendix we calculate the number n(k) of planar graphs with k internal $-lines. 
This number appears in the expression ( 3.2 1) ) for the gluon self-energy. As we already 



mentioned in section 3.2, in a planar graph both ends of a <l>-line attach to the gluon line 
either from above or from below. Let there be a loops above and b loops below the gluon 
line, so that 

k = a + b (B.l) 

The number C a of ways to combine the 2a vertices with a planar loops above the gluon is 
independent of b, and is the same as for the quark propagator. If each loop gives the same 
weight x to the quark propagator S(x) (ignoring its Dirac structure) then 



S{x)=Y,C a x a (B.2) 



a=0 



The function S{x) is determined by the DS equation which generates all planar diagrams 
for the quark propagator (c/. Fig. ||), namely S(x) = 1 + xS 2 (x) with S(0) = 1. This gives 



S(x) = (B.3) 

A Taylor expansion shows that S(x) is the generating function for Catalan numbers. Thus 
the numbers C a in ( |B.2[ ) are 



a + 1 \ a 

For a given ordering of the 2a + 2b vertices on the gluon line we have then C a Cb ways 
of drawing planar loops. The number of distinct orderings of the vertices is given by the 
number of ways to choose 2a vertices (regardless to their order) from a set of 2k vertices, 
i.e., by the binomial factor (^). The total number n(k) of distinct planar diagrams with 



k loops to be used in ( 3.21| ) is thus: 

k /2k\ 

n(k)=J2L )C a C k -a = C k C k+1 (B.5) 

a=0 ^ a ^ 

where the last equality may be derived as follows. Rearranging factors in ( |B.5[ ), 

a=0 1 \ / \ / a=0 V / V 1 / 
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Comparing the coefficients of x + in the equivalent expressions 

n=0 a=0 v 7 v 7 

(B.7) 



we get 




(B.8) 



Using this in ( |B.6| ) we obtain 

C. Decoupling of zero-momentum gluons from color-singlet quark loops 

In section [O] we gave a formal proof that any number of external zero-momentum pho- 
tons decouple from an electron loop contribution to the photon self-energy, due to charge 
coherence. Here we shall extend this proof to the decoupling of zero-momentum gluons 
from a (color singlet) quark loop correction to the photon propagator. 

Let n^f... Mn (p) denote a quark loop correction to a pho- 
ton propagator with momentum p and Lorentz indices a, /3. 
The n external, zero-momentum gluons attached to the quark 
loop have Lorentz indices \x\ . . . (J, n , and their color indices 
a\, . . . ,a n are implicit. The QED result 

n^»=0 for n>l (C.l) 

also holds in QCD for n = 1 since TrT Ql = 0, and for n = 2 
because the color factor is then simply oc 5 aia2 . For n > 3 
the QCD amplitude LT^f... Mn (p) differs from the QED one 
because the various diagrams are weighted by different color 
factors of the type Tr [T ai . . . T a ™]. We now show that the 
part of the QCD amplitude which corresponds to a given 
color factor formally vanishes by itself. This ensures the 
vanishing of the complete QCD amplitude. 

Consider the part of the corresponding QED amplitude 
which is built from all diagrams that have the same cyclic 
ordering of the external photon lines [i\ . . . fj, n , where the indices are numbered by following 
the fermion loop against the direction of the fermion arrow. There are many diagrams of 
this type, which are distinguished by the positions of the two photon lines with momentum 
p and indices a, j3 in the sequence of cyclically ordered external photons (see Fig. We 




Figure 7: A QED diagram 
where the photons with van- 
ishing four- moment a (denoted 
by dashed lines) have a given 
cyclic ordering fii . . . fi n . Lines 
Ui and u n are in the proof as- 
signed momenta I ^ 0. 
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temporarily take the photons with indices fi± and \i n to carry incoming momenta / and 
—I, respectively, while the photons (12 ■ ■ ■ fJ* n — 1 have vanishing four-momenta. Denoting the 
amplitude just described by ^^..^(p, I), we wish to show that T^f... lln {p,0) = 0. 

The momentum I may be taken to flow in the direction of the fermion arrow between the 
vertices Hi and /i n , and not beyond. According to ( ft.3| ) the derivative —ed / dl^ n+1 applied 
to T^,,^ n (p,l) then generates all QED diagrams with (n + 1) external zero-momentum 
photons having the cyclic ordering \x\ . . . fi n +i- 



- T alS (v I) 

- L fi 1 ...fi n \P'i') 



Assuming that the Taylor expansion 



= ^f..^ n+1 (p,Z = 0) (C.2) 

1=0 



+ 0{l 2 ) (C.3) 

1=0 



is well-defined we obtain 

Vt., n (P, = Vl, n (p, 0) - \ T^ n+1 (p, 0) + O (/ 2 ) (C.4) 
It is readily seen from the charge conjugation symmetry that 

22f..^ 1 (p.0 = (-i) B 2g?.^.(p,0 (c.5) 



Reversing the order of //1 . . . /i n in ( |C.4| ) and recalling that the amplitudes were constructed 
to be cyclically symmetric in the Lorentz indices of the zero-momentum photons, 

T£..^ n+l (P, 0) = T«f +lMn ... Ml (p, 0) (C.6) 

we get 

(-i)"2g!.^(p.O = (-i)%f..M»(^°) - Vir +1 ^r;^„ n+i feo)+o(! 2 ) (c.7) 



Comparing ( |C.4| ) and ( |C.7| ) we arrive at 

2g?.^ fl (p,0) = (C.8) 
for n > 2. Since all QCD amplitudes with the same cyclic ordering \i\ . . . \i n of the external 



gluons have the same color factor Tr [T ai ...T an ] this proof for QED implies that (C.l) 
holds also in QCD. 

D. Quark-photon vertex T^(k,k) 

In this appendix we solve the implicit equation ( |4.7| ) for the 7qq vertex in the case of 
the chirally invariant quark propagator (|3.7|), 



Si(p) = a p p ; a p = ^ ll - J 1 - (D.l) 
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Using this expression for S(p) ( [4.7] ) becomes 

where we introduced the dimensionful parameter /, 

/ = l?a k a- k (D.3) 
Chiral and parity invariance restricts T^(k,k) to the form 

T"(k, k) = A)7 M + Mk^ + A 2 k^ + A 3 k"^t + A 4 k^ + ^5756^(7, k, k) (D.4) 
where we defined 

e»( 1 ,k,k) = e^ lv k p k CT (°- 5 ) 
We find the coefficients Ai by inserting flD.4j ) into ( |D.2| ) and using 

= k^ + k^ty-k-k^ - ij 5 e^(j, k, k) 
^75^(7, fc, % = -i7 5 e^( 7 , k,k)k-k + y> 1 [k 2 k 2 - (k ■ k) 2 ] 

-kv[k 2 k- -k-klf:]-k^[k 2 ^ -k-k\f\ (D.6) 



This gives the conditions 

Ax = fk 2 (A 2 -A 5 



A = 1 - fk • k {A + k ■ k A 5 ) + fk 2 k 2 A 5 

2/ 



with solutions 



A 2 = KAo + k-kAs + ^Ax) 
A3 = f{A Q + k-kA h + k 2 A A ) 
A A = fk 2 (A 3 -A 5 ) 

A 5 = -f(A + k-kA 5 ) (D.7) 



1 + 2fk ■ k + f 2 k 2 k 2 ' 1 + 2fk ■ k + f 2 k 2 k 2 

i2a - 2a 2fk 2 k 2 A A A l + f 2 k 2 k 2 A 

k 2 A 1 = k 2 A, = - i J l f2k2 - k2 A 5 ; A 2 = A 3 = - x _ J pk2 - 2 A 5 (D. 

The result for T fl (k,k) then follows from ( p.4| ): 



1 + 2fk ■ k + / 2 £: 2 A; 2 

f2 f/1 , f2i2l.2\ 



+ i- 2 /w (fc ^ 2 + ^^ 2) + f i-pv& + (D ' 9) 

Let us check that this expression for the vertex satisfies the Ward-Takahashi relation 
[Q|), Straightforward algebra yields: 



1 _ ft-2 1 _ fl-2. 

^ ( ">=*T374*-*w4* (D - io) 
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Using flOp and Eq. Q (for 6 = 0), 



2 2 

fi a p 



(D.ll) 



we get 



i-A 2 



l-f 2 k 2 k 2 k 2 a k ' l-f 2 k 2 k 2 k 2 c 



(D.12) 



Substituting these expressions in (D.10) gives the Ward-Takahashi relation ([4.9| 
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